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1. Introduction 

Potential models of hadrons originate from the quantum chromody- 
namics (QCD), but they are not rigorously deducible from the the- 
ory. Rather, these models are substantiated by various approximate 
approaches and estimates in QCD.0 Various models have their own 
areas of application. In particular, the linear potential which follows 
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from area law in the lattice approximation of QCD, describes, by 
construction, the static interaction of quarksB Thus it can lawfully 
be exploited only in nonrelativistic potential models J 

The description of light meson spectroscopy needs the develop- 
ment of appropriate relativistic models. They frequently are built 
as single-particle wave equations@ _ which is not satisfactory. Actu- 
ally, mesons should be treated as composite two-quark relativistic 
systems. The reliable basis for this purpose is the relativistic di- 
rect interaction theory (RDIT) presented by various approaches and 
formalisms,^ such as Lagrangian formalism with higher derivatives,^ 
relativistic Hamiltonian mechanics canonical formalism with 
constraints,^ Fokker-type action formalism^ etc. 

Given a nonrelativistic potential, RDIT determines the general 
structure of its relativistic counterpart. In so going the great arbi- 
trariness arises in the choice of concrete relativistic interaction. Con- 
sequently, the variety of relativistic potential models has appeared in 
the literature.IlM! Each of them has its own advantages and areas of 
application, but these models are not substantiated by QCD better 
than nonrelativistic models. 

A possible way to substantiate relativistic direct interactions leads 
through classical field theory. Especially, we mean the Fokker-type 
action formalism which, among other approaches to RDIT, is most 
closely related to this theory. There exists the class of Fokker actions 
which correspond to particle interactions via linear fields, such as 
scalar, vector, and other tensor fields.^ These actions are built on 
the solutions to relevant wave equations. In the nonrelativistic limit 
they lead to the same Coulomb (or Yukawa) potential. 

Working within this scope for the confinement case, one could 
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But no such solutions leading to confining potentials are known in 
the literature. Moreover, they are believed to not exist due to the 
essentially quantum nature of confinement. This is concerned with 
standard Yang-Mills theory while there exist various nonstandard 
theories which involve effective Yang-Mills fields arising from QCD. 
These theories may be used as sources of confining potentials. 

In the present paper we find the relation between certain non- 
standard classical theory of Yang-Mills field and the Fokker-type 
confinement model. The former is developed in Ref. |T7|. This theory 
describes non-Abelian gauge field averaged over quantum fluctua- 
tions. It is based on the effective Lagrangian obtained from the 
study of infrared behaviour of gluon Green's functions in QCD.Etl 
Field equations following from this Lagrangian are of 4th order, and 
some static non-Abelian solutions to them have been used in a sort 
of bag confinement model.^ 

Here we obtain from this theory the Abelian retarded and time- 
symmetric potentials of moving pointlike source. Both of them are 
of confining type and reduce in the nonrelativistic limit to the lin- 
ear potential. Then, using these potentials, we construct the time- 
asymmetric and time-symmetric Fokker-type actions. The latter is 
already known in the literature. Two equivalent versions of this ac- 
tion have been proposed by Rivacobal^ and WeissS It is noteworthy 
that both the authors proceeded from general preliminaries of RDIT, 
without referring to field-theoretical interpretation of particle inter- 
action. 

The time-symmetric action leads to difference-differential equa- 
tions of motion which are difficult to deal with. The only circular- 
orbit solutions to these equations are found in Refs. [2D] and ^T]. Con- 
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is well defined in terms of second-order differential equations of mo- 
tion. Thus this action can be considered as the classical background 
model of relativistic two-quark quantum dynamics. Following^ we 
reformulate this model into the Lagrangian formalism. Then we 
transit to the Hamiltonian formalism, and integrate the two-body 
problem in quadratures. 

The time-asymmetric analogue of Rivacoba- Weiss model is the 
simplest version of relativistic confinement model. It can be appro- 
priate for the classical description of light mesons for which the con- 
finement interaction dominates. To include into consideration also 
heavy mesons one can modify the present model by adding to the 
action the vector-type interaction term from the time-asymmetric 
version of the Wheeler-Feynman electrodynamics. This corresponds 
to the taking account of Abelian solution to the standard Yang- 
Mills equations (i.e., the classical analogue of one-loop correction 
in QCD). In the nonrelativistic limit this mixture leads to the well 
known Coulomb plus linear potential. The modified model becomes 
appreciably cumbersome but still remains solvable. 

Here we do not propose a quantum version of the present model. 
Instead, we make some estimates of the Regge trajectory from clas- 
sical and semiclassical considerations and obtain a physically reason- 
able result. 

The paper is organized as follows. In Section 2 we obtain the 
Abelian potentials of moving pointlike sources from the standard 
and nonstandard theories of Yang-Mills field. The formers are the 
Lienard-Wiechert potentials and their causal modifications while the 
latters turn out to be the modifications of potentials proposed by 
Weiss. They are obtained with the Green's functions found in Ap- 
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ing from the standard and nonstandard theories, and construct cor- 
responding time-symmetric and time-asymmetric Fokker-type inte- 
grals. The latter is used as the base of time-asymmetric confinement 
model. In Section 4 this model is reformulated in the framework 
of the Hamiltonian formalism. Various special cases of two-body 
problem are considered in Subsections 4.1-4.3 and Appendix B. Es- 
timates of Regge trajectory are quoted in Section 5. Section 6 is 
devoted to general discussion of the model. 

2. Abelian potentials from the standard and 
nonstandard theories of Yang-Mills field 

We shall consider both the standard and nonstandard classical the- 
ories of the Yang-Mills field. The standard theory (ST) is based on 
the well known Yang-Mills Lagrangiani 



The nonstandard theory (NT) proceeds from the effective 



Here the components of the gauge field /K fl (x) (/i = 0,3) and the 
current of sources J /i (x) take values in the Lie algebra Q of gauge 
group; [X, Y], (X, Y), and V M X = d^X— [A M ,X] are the Lie brackets, the 
Killing-Cartan metrics, and the covariant derivative, respectively, 
defined for any X, Y G Q\ = d^k v — d v k^ — [A M ,Aj,] is the tension 
tensor; k, is some parameter of the dimension of inverse length, and £ 
is meant here as an arbitrarv dimensionless narameter. We snnnose 
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Killing-Cartan metrics is nondegenerate and positively denned. In 
the real matrix representation it can be presented in the form 

(X,Y>=-— tr(XY), 

where the number Nr depends on the representation chosen. Greece 
indices move due to the metrics 77^ of the Minkowski space-time IM4 
which is chosen timelike, i.e., \\rj^\\ = diag(+, —,—,—). 

Field equations following from the Lagrangians (1) and (2) are 

V„F^ = 4ttJ^, (3) 

and 

{2V,V 2 - (1 + VaV,V a + <eV 2 V,} F^ 4 = 4ttk;V, (4) 

respectively. Both of them are compatible provided the current ^ is 
covariantly conserved, 

= 0. (5) 

In the present paper we are interested in the relativistic system 
of N pointlike charged particles interacting via the Yang-Mills field. 
The current J M corresponding to this system is^ 

= E J2C&) = E fdr a Q a ^5(x - z a ). (6) 

a a J 

Here z%(r a ) (fi = 0,3, a = 1,N) are the space-time coordinates of 
ath particle world line in IM4 parametrized by an arbitrary evolu- 
tion parameter r a , z%(r a ) — dz^/dr a , and Q a (T a ) is the charge of ath 
particle. Substituting (6) into (5) one obtains the Wong equations^! 
determining the evolution of charges, 

Qo = i?[A AI (« ),Qa] J a = TJT. (7) 

The total action corresponding to field + particle system can be 
written down as follows: 
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where C is £st or £nt, and m a is the rest mass of ath particle. 
The variation of the action (8) over A M yields the field equations 
(3) or (4). Varying this action with respect to particle positions z% 
and taking account of (7) one can obtain the following equations of 
particle motion: 



In order to determine motion of particles it is necessary to solve 
the total set of linked equations, namely, the field equations [(3) or 
(4)], the Wong equations (7), and the equations of motion (9). 

We intend to formulate a particle dynamics in the scope of RDIT. 
For this purpose one should eliminate field variables A^x) in favour 
of their expressions in terms of particle positions z% and, possibly, 
charges Q a . In other words, it is necessary to find a solution to field 
equations. But this task is very complicated because of nonlinearity 
of the problem. 

Here we limit ourselves by search of Abelian solutions to field 
equations. Let us suppose that 



A„{x) = nA^x), f{x) = n J"(x), Q a (r a ) = nQ a (r a ) (10) 



etc., where n is a unit constant vector in Q. In this case all Lie- 
bracketed expressions vanish, in particular, 

F^(x) = nF^(x), F^(x) ee d^A v (x) - dvApix), (11) 

and V M reduces to d^. The Wong equations (7) yield 




d mz, 




(9) 



Qa = 



Qa = q, 



a 



= const. 



(12) 



Then the field equations (3) and (4) reduce to 
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and 

d v UF^ = 4ttk 2 J^ (14) 
respectively, where the current 

J»(x) = £ J£(z) =T,q a [dr a z^(x - z a ) (15) 

a a J 

is conserved identically, i.e., <9 M J M = 0. Due to this fact both the 
equations (13) and (14) are gauge invariant with respect to the one- 
parametric (compact) group of residual symmetry. 

At this point we have come to the linear field equations which can 
be solved by means of the Green's function method. In the standard 
case we deal exactly with the electromagnetic problem. Using the 
Lorentz gauge fixing condition, 

d^A" = 0, (16) 

we reduce the equation (13) to d'Alembert equation, 

DA" = 4ttJ^, (17) 

and immediately obtain its solution, 

A» = D ri * J", (18) 

where * denotes the convolution, and 

D v (x) = (l + r]sgnx )5(x 2 ), (19) 

is one of the retarded (rj = +1), advanced (rj = —1), or time- 
symmetric (77 = 0) Green's functions of d'Alembert equation. 

Let us consider the equation (14) of the nonstandard theory. Us- 
ing the Lorentz condition (16) one reduces it to the following equa- 
tion: 
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which is of 4th order. In Appendix A the corresponding retarded, ad- 
vanced, and time-symmetric Green's functions are calculated. They 
are: 

E v (x) = lK 2 (l + r]sgnx )G(x 2 ). (21) 

Thus the solution to (20) reads as (18), but with E n instead of D v . 

Actually, the linearity of equations (17) and (20) allows solutions 
of more general structure, 

^ = E^ = £GW£, (22) 

a a 

where G Va = D Va for ST, and G Va = E J]a for NT. Here rj a take values 
+1, -1, or 0, each own for different particles. 

In an explicit form the solutions (22) can be written down as 
follows: 

A"(x) = £4f(*) =T,qaJdrXG Va (x - z a ), (23) 

a a J 

where the quantity A%(x) represents the relativistic potential created 
by ath particle. In both the ST- and NT-cases each particle poten- 
tial (as well as the total sum (23)) satisfies the Lorentz condition 
(16). 

Up to the numerical factor, the only difference between (19) and 
(21) is that the function 5(x 2 ) is replaced by Q(x 2 ). This substitution 
was guessed by Weiss in Ref. £T] where the time-symmetric potential 
(in our case, Eqs (23) with G Va = Eq, a = 1,N) has been proposed 
for the model of the action-at-a-distance linear confinement. 

3. Equations of motion and Fokker-type action 
integrals 

Now the equations of particle motion can be obtained in a closed 
form bv substitution of the relativistic ootentials and the con- 
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standard case this procedure leads to an appearance of divergent 
self-action terms which can be regularized in usual way.S The re- 
sulting equations of motion can be presented in the form: 

± % = q a E F^zl + R afl , (24) 

where 

Fabvv = 2Qb Jdr b (l + rj b sgn z Q ab )8' (z 2 ab ) {z ab[L z bv - z abv z b ^} , (25) 
z a b = z a - z b , and 

" 3™' 1 "~^T| *^ *^ ' <26) 



The self-action terms i£ a/i correspond to radiation reaction. They 
disappear if fields generated by particles are time-symmetric (i.e., if 

Va = 0). 

In the nonstandard case no divergences and self-action terms arise. 
Thus the equations of motion are calculated immediately. They are 
described by (24) with R afi = and 

Fabiiv = \^qb Jdr b (l + T] b sgn zl b )8{z 2 ab ) {z ab ^z bv - z abiy z bfi } . (27) 

We have obtained the closed set of equations of particle motion 
which are not obvious to be directly deducible from the variacion 
principle. Below we construct the relevant Fokker-type version of 
the theory and examine it consistency with the equations obtained 
above. 

The purpose is to eliminate field variables from the total action 
(8). Using (10)-(12) in (1), (2), and then in (8), one obtains the 
action 

I = -ffree + -fiat + Afield 5 (28) 

where 
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lint = ~fd A x J»A„ (30) 
are the same for ST and NT while /fi e id is different: 

Afield = 

-—Jd 4 xF^F^ (31) 



16tt* 
for ST, and 

/field = Jd 4 x (d x F^)(d x F^), (32) 

for NT. The term ig c id can be transformed to the form 

SF/^W+fSD' (33) 

where H"" = F v » for ST, and H"" = DF^/k 2 for NT. Taking into 
account the field equations (13) and (14) in r.h.s. of (33) and omit- 
ting surface terms, we obtain 

I = /free + |/ int . (34) 

Now substituting the current (15) and the potential (23) into (30), 
one can present the second term in r.h.s. of (34) in the following 
form: 

;Kint = E E + I E ^aa, (35) 

a < b a 

where 

hb = ~q a qb J [dr a dr b z a ■ z b G Vba (z ab ) , (36) 

and r\ ba = \{r) b — rj a ). In the ST-case the self-action term I aa di- 
verges. It can be regularized and unified with ath term of /f ree - In 
the NT-case this term vanishes. Thus in the both cases the resulting 
interaction term |ij nt has the form: 

\hnt = - E E Wb f fdr a dT b z a ■ z b G Vba (z ab ) , (37) 

a < b 

where G Vba = D Vba for ST, and G Vba = E Vba for NT. 
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In the case of NT each constituent (36) by means of integration 
via parts (see Ref. [25] for such a technique) can be transformed to 
the following form (here we omit all unessential constant factors): 



-oo +oo 



J I dr a dT b z a -z b (l-\- T] ba sgn z Q ab )Q(z 2 ab ) 

-oo — oo 

+oo +oo 

= -2 J J dr a dr h (z ab ■ z a )(z ab ■ z b )(l + rj ba sgn z° ab )5(z 2 ab ) 

o 

(38) 



oo — oo 

-| r a =+oo-| T fc = 

(l + 7ft a sgnz° b )Q(zl b )z 2 ab 



-oo. 



7V 



The second term in r.h.s. of (38) is divergent, but it does not con- 
tribute in equations of motion and can be omitted. Then the inter- 
action term (37) for NT can be put in the equivalent form, 

|4it = ~7T EE Qaqb I dr a dr b (z ab ■ z a ) (z ab ■ z b )D Vba (z ab ) . (39) 

^ a < b 

Fokker-type equations of motion following from this action dif- 
fer from those (24) directly obtained from the field theory. Firstly, 
they do not reproduce the self-action terms R a/I which, in general, 
are present in r.h.s. of equations (24) for ST. In this paper we sup- 
pose that these terms can be neglected since in QCD a radiation 
is suppressed by confinement. Secondly, the sign factors j] b in the 
expressions (25) and (27) for F ab[lv are replaced by rj ba . This changes 
the causal structure of pair particle interactions. Namely, while equa- 
tions (24) correspond to retarded, advanced, or time-symmetric fields 
generated by 6th particles (and sensed by ath particle) for j] b = +1, 
-1, or 0, respectively, in the Fokker-type equations the causality of 
interactions is its own for different pairs of particles. 

There are only two cases in which the direct interaction can be 
treated as a field-type one. The first case corresponds to the time- 
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the action (34), (29), (37) in this case coincides with the Wheeler- 
Feynman action of time-symmetric electrodynamics.^ For NT it cor- 
responds to the action-at-a distance confinement model in the form 



The second case which is tractable in terms of field interaction 
realizes only for two-particle systems. It corresponds to the choice 
rj2 = —rji = 7721 = 77 = ±1. For ST this is the case of the time- 
asymmetric electromagnetic interaction proposed by Staruszkiewicz, 
Rudd and Hill,^ and studied in more detail in Ref. |28|. For NT the 
corresponding time-asymmetric Fokker-type action can be taken as 
the classical base for relativistic confinement model. 

4. Time— asymmetric model with confinement 
interaction 

In this section we consider the two-particle model available for the 
classical description of mesons. It is based on the time-asymmetric 
Fokker-type action which combines interaction terms (37) from ST 
and (39) from NT. Since mesons are chargeless systems, we put q\ = 
— (?2 = q- Then the time-asymmetric Fokker-type action has the 
form: 



where a = q 2 , (3 = ^q 2 K 2 , and r\ = ±1. In the nonrelativistic limit this 
action leads to the well known interquark potential U = —a/r + (3r. 
Integrating the second term of the action (4) once, we reduce 




/ = 




x {az r z 2 - (3(zi2 ■ £i)(zi2 • z 2 )} , 



(40) 
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of our model in the framework of a manifestly covariant Lagrangian 
formalism with the Lagrangian function 

L = 6F(a U (j 2 ,S), (41) 

where = rjy-z > 0, z = Zi - z 2 , y = (z\ + z 2 )/2, a a = \fz~l/0 > 0, 
5 = Z\ • z 2 /0 2 > 0, and with the holonomic constraint z 2 = 0, rjz° > 
0. All variables in (41) depend on an arbitrary common evolution 
parameter r. In our case the function F has the form: 

F = £ m a a a -a5 + f3. (42) 

a=l 

We note that quantities 9,a a ,S in r.h.s. of (41) and (42) are well 
defined and positive if particle world lines are timelike. 

The transition to the manifestly covariant Hamiltonian descrip- 
tion with constraints leads to the mass-shell constraint which deter- 
mines the dynamics of the model and has the following form:^ 

0(P 2 , V 2 , P-Z, VZ)= fre e + </>int = 0. (43) 

Here = — z^P -w / P • z; P^ and are canonical momenta 
conjugated to y^ and z^, respectively; the function 

free = \P 2 - \{m\ + m 2 ) + (ml - + v 2 (44) 

corresponds to the free-particle system, 



a(P 2 — ml — m 2 ,) a 2 ml 
0int = „ 1 n — 



rjP-z rjP-z a =ib a + a 

describes the interaction, and 

b a = r)(\P-z+ (-fv-z) , a =1,2, a = 3- a. (46) 

We note, that the quantities a a are related to canonical variables 
by the equations: 
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Since a a must be positive, the following conditions arise: 

6 fl + a>0, a =1,2. (48) 

They restrict the whole phase space to a physical domain in which 
the Hamiltonian description is equivalent to the Lagrangian one. 

In order to study the dynamics of the present model it is conve- 
nient, following Ref. |22], to transit from the manifestly covariant to 



three-dimensional Hamiltonian description in the framework of the 
Bakamjian-Thomas model Within this description ten genera- 
tors of the Poincare group P M , as well as the covariant particle 
positions z% are the functions of canonical variables Q, P, r, k. The 
only arbitrary function appearing in expressions for canonical gener- 
ators is the total mass \P\ = M(r, k) of the system which determines 
its internal dynamics. For the time-asymmetric models this function 
is defined by the mass-shell equation^ which can be derived from 
the mass-shell constraint via the following substitution of arguments 
on the l.-h.s. of (43): 

P 2 ^M 2 , <u 2 ^-k 2 , P-z^r]Mr, v-z^-k-r; (49) 

here r = |r|. 

Due to the Poincare-invariance of the description it is sufficient to 
choose the centre-of-mass (CM) reference frame in which P = 0, Q = 
0. Accordingly, Pq = M, Jqi = (i = 1,2,3), and the components 
Si — \si Jjk form a 3-vector of the total spin of the system (internal 
angular momentum) S = r x k which is the integral of motion. At 
this point the problem is reduced to the rotation invariant problem 
of some effective single particle; such a problem is integrable in terms 
of polar coordinates, 
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Here S — |S|; the unit vectors e r , are orthogonal to S, they form 
together with S a right-oriented triplet and can be decomposed in 
terms of Cartesian unit vectors i, j: 

e r = icos<^ + j sin if, = — isin^ + j cos<^, (51) 

where ip is the polar angle. 

The corresponding quadratures read: 

t - t Q = Jdr dk r (r, M, S) /dM, (52) 
ip - (f Q = -Jdr dk r (r, M, S) /dS, (53) 

where t = \{z® + z^cm is the fixed evolution parameter (unlike the 
undetermined one r), and the radial momentum k r , being the func- 
tion of r, M, S, is defined by the mass-shell equation written down 
in terms of these variables, 

(J)(M 2 , -k 2 , rjMr, -k • r) 
- |m 2 , -k 2 r - ^, 7?Mr, -/c r rj = 0. (54) 

The solution of the problem given in terms of canonical variables 
enables to obtain particle world lines in the Minkowski space using 
the following formulae:^ 

s2 = * + K-)V, (55) 

k fl - k \ S 
z a = K-)° r + + TGr + ^M e ^- ( 56 ) 

Especially, the vector z = zi — Z2 = r characterizes the relative 

motion of particles. 



4..1 Purely confinement model 

Hereafter we restrict ourselves to the svstem of eaual rest masses. 
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interaction. The mass-shell equation in this case reads: 



It easy to obtain from (57) the expression for k r (r, M, S) 




(57) 



k r = eJf(r,M,S), e = ±1, 



(58) 




(59) 



Besides, we must take into account the condition: 



\M 2 -k 2 r >0 



(60) 



which follows from (48). Then from (58)-(60) we obtain the restric- 
tion: 



The quadratures (52), (53) with (58), (59) can be reduced to the 
elliptic integrals. Here we omit their expressions. The integration 
is spread over the domain of possible motions (DPM) which is de- 
termined by the conditions (59) and (61). In the case S > DPM 
consists of the connected interval r\ < r < r2, where ri, are posi- 
tive roots of the equation f(r, M,S) =0. The latter can be presented 
as the reduced cubic equation with respect to 1/r: 



It has two real positive solutions provided the following conditi 
holds: 



< r < \M/(5. 



(61) 




(62) 



M > M C {S) 



(63) 



where the function M C (S) is defined in the implicit form 



/\/f2 / 



/l^2\ 3/2 
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The equality in (63) corresponds to the case T\ = = r c of circular 
particle orbits with the distance between particles 



satisfying the set of equations: 

/(r c , M c , S) = 0, 0/(r„ M c , S)/dr c = 0. (66) 

In the limit S — > the quadrature (53) yields ip = ipo, and a par- 
ticle motion becomes one-dimensional (i.e., in the two-dimensional 
space-time IM2 parametrized with x° and, say, x 1 ). Besides, T\ — > 0. 
Thus DPM becomes < r < r<i- The point r = corresponds to 
particle collision. This point is not singular for the quadrature (52) 
and particle coordinates (55), (56). Thus the motion of particles can 
be smoothly continued as if they pass through one another. 

4. .2 General model, S > 

Let us consider the general case a > 0, (3 > 0. The corresponding 
mass-shell equation can be written down as follows: 



A ( a a 2 \ f / 3r\ A m 2 (M a 
-jA - I A - -jj J M - 2^:j A - 2— f — + — 



It is quadratic equation with respect to 

2 



/M ay 9 , x 

A=( T + -J -^ 2 >0, (68) 

where A must be positive because of the conditions (48). As to k r 
the equation (67) is biquadratic. Its solution can be presented in the 
following form: 

k r = e^Jf ± (r, M, S), e = ±1, (69) 
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where 

f ±i r,M, S) - ^ + M^-^§, (70) 

/i±(r,M,5) = p(r, M, S) T \jd{r, M, S) (71) 

d(r, M, 5) = £ 2 (r, M, 5) + (l - 2^) (72) 

, ™ ™ 2 ( M a \ f, a 2 S 2 

ff (r, M, 5) = - ( T + 7 ) - (l - 2^) ^ + ^. (73) 

Among two solutions /± for k 2 we choose that one which is smooth 
in DPM and reduces to / (59) in the limit a — > 0. 

DPM is analyzed in Appendix B. In the case S > we have 
d(r,M,S) > 0, r > 0. Thus both the functions f±(r,M,S) are 
smooth provided r 7^ ±M/(3, and f-(r,M,S) reduces to / (59) in 
the limit a — > 0. DPM in this case is determined by inequality 
f-(r,M,S) > provided the condition (68) holds. Similarly to the 
purely confinement case we obtain r\ < r < where r\, are 
positive roots of the equation: 

S 2 (M a\ / ~ «\ 

X {{ 1 - 2 M){Y + r)- 2 M\ = (74) 
It can be reduced to a fourth-order algebraic equation which has 
two real positive solutions provided (63), where the function M C (S) 
can be presented in a parametric form, 

M?(X) = 2ma A[4 + 5A + 2A' + y (l + A)»(4 + 3A)] 



(1 + A) 
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52(A) = a2 A(2 + A)»[l + „(1 + A)*] 



4 + 5A + 2A 2 + z/(l + A) 2 (4 + 3A) ' 

v = — ± < A < 00. 75 

m z 

The condition M = MJS) corresnonds to circular orbits with 
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Figure 1: General model, S > 0. Typical example of phase (left) and particle (right) 
trajectories. Arrows show the direction of evolution; H — centre of mass. 



Our attempts to express the quadratures (52), (53) with (69)-(73) 
in terms of known (elementary and special) functions have not met 
with success. Thus we calculated them with a computer. Neverthe- 
less by means of analytic calculations it can be shown that particle 
world lines in IM4 are timelike and smooth curves. They represent a 
bound motion of particles for all values of M allowed by (63), (75), 
and S > 0. The typical example of phase and particle trajectories 
are shown in figure 1. 

4. .3 General model, S = 

In the case S = we have d(r, M, 0) = g 2 (r, M, 0). Since there 
exists the point ro < a/m such that g(ro,M, 0) = 0, the functions 
/±(r, M, 0) are not smooth. Moreover, in the domain r < tq the 
function /_(r, M, 0) has not the proper form in the limit a — > 0. 
Thus the evolution of particles cannot straightforwardly be continued 
farther. 

We point out that the distance ro at which the smoothness of 
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of the classical electron radius. In the case of strong interaction the 
distance ro and the Compton length of quarks can be commensu- 
rable quantities. Thus the classical description of particle motion 
at r < ro may be important for the construction of quantum the- 
ory. Especially, this is concerned with the case of S-states. Below 
we propose the way to continue the particle motion in the domain 
r < ro- It leads beyond the rigorous treatment of analytical mechan- 
ics and therefore cannot be a reliable basis of quantum-mechanical 
description. But it will be noted that this method arises naturally 
from the present model itself. 

Let us choose in r.h.s. of (69) the function: 

/o(r,M) = I 

l/_(r,M,0), r > r 

(M a\[M a 2m 2 /M } , . 

= {-2 + ^){-2 + -r-T^WlMl <76) 
which is smooth provided r ^ \Mj (3, and reduces to (59) if a — > 0. 
DPM in this case is < r < while the point r = is critical: 
A — > 00, r — > 0. This means that the equivalence between the 
Lagrangian and Hamiltonian formalisms violates. It can be shown 
that at the collision one of particles reaches (but not exceeds) the 
speed of light while another does not. Again, the particle world lines 
should be somehow continued farther. 

The existence of a singular collision point is due to the time- 
asymmetric vector interaction from ST. The confinement interac- 
tion does not change qualitatively the behaviour of particles in the 
neighbourhood of collision. Specific features of the time-asymmetric 
model with attractive vector interaction in IM2 have been analyzed 
in Ref. |^T]. Following this work, in the framework of Hamiltonian de- 
scrintion the interest,! ner nossibilitv exists to continue smoothlv world 
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Figure 2: General model, S = 0. Typical example of phase trajectory continued in 
the non-Lagrangian domain r < 0. Arrows show the direction of evolution. 




Figure 3: General model, S = 0. Typical example of particle world lines. LS - 
Laexaneian segments: NS - non-Laeraneian segments: • - collision: one of Darticles 
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purpose to suppose formally that after the collision the variable r 
becomes negative. Then we have /o(r, M) > 0, r G [— 2a/M, 0), 
and the motion of particles can be continued up to the distance 
|r| = 2a/ M. At this point which is also critical both particles reach 
(but not exceed) the speed of light. Again, one can smoothly con- 
tinue world lines up to the next collision e.t.c. We note that A < 0, 
r G (— 2a/ M, 0). Thus the segments of world lines obtained as above 
do not follow from the Lagrangian description. The resulting world 
lines combine the Lagrangian and non-Lagrangian segments sepa- 
rated by the collision points. They describe the bound periodic mo- 
tion of particles. The corresponding phase trajectory and world lines 
are shown in figure 2 and 3 respectively. 

The formal continuation of evolution proposed above permits 
some reinterpretation in terms of the Lagrangian description. Ex- 
pressing the quantities 0, a a and 8 in terms of canonical variables 
one can examine that some of them have wrong (i.e., negative) sign 
if r < 0, i.e., if particles pass non-Lagrangian segments of world lines. 
Equivalently, one can keep r > changing signs of some constants 
m a , a, and (3 in the Lagrangian (41), (42). In such a manner one 
can realize that particles move as if each one changes signs of its rest 
mass and charge, m a — > — m a , q a — > — q a , once it passes a critical 
point with the speed of light. We note that at the collision point one 
of world lines is timelike. Thus the mass and charge of this particle 
remain unchanged up to the next critical point. As a result, after 
having passed the non-Lagrangian segment each particle returns its 
proper values of mass and charge. 
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5. Semiclassical estimates of Regge trajectory 

It is well known that nonrelativistic potential model with the lin- 
ear potential leads to the Regge trajectory with the unsatisfactory 
asymptote M ~ S 2 ^. Here we do not propose a quantum version of 
the present model, but we make the estimates of the Regge trajectory 
from what follows. 

Usually the Regge trajectories in the potential models are cal- 
culated in the oscillator approximation^ Then the leading Regge 
trajectory originates from the classical mechanics: it is close to the 
curve of circular motions on the (M 2 , S^-plane. 

In our case this curve follows from the equation (64) for a = or 
from (75) in the general case. The latter has in the ultrarelativistic 
limit M c — > oo the desirable linear asymptote: 

M 2 « 6VS(3S + 6 (ra 2 - Sa/3) . (77) 

It is remarkable that this asymptote is achieved only by taking ac- 
count of a relativity. 

Let us compare the classical Regge trajectory (64) of purely con- 
finement time-asymmetric model to that which follows from the time- 
symmetric Fokker-type confinement model with the same parameters 
mi = 77i2 = m and (3. Considering the circular orbit solution, given 
in Ref. |20|,B for large M one can obtain: 

M 2 w 4(1 + sintf){2costf/3S + m 2 }, (78) 

where the angle 

•d w 0.7391 « 0.2353tt (79) 
is the solution of the transcendental equation 

tf = costf. (80) 
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(3S 10 




Figure 4: Classical Regge trajectories; - - from general time-asymmetric model 
at various rates of parameters; - - from time-symmetric model. 

The only difference between the asymptotes (77) (with a = 0, of 
course) and (78) consists in slightly different numerical factors at 
the linear and constant terms 

linear terms constant terms 

6^ « 10.3923 6 (81) 

8 cos (1 + sin ) « 9.8955 4(1 + sin ) « 6.6944 

Moreover, it turns out that by the substitution # = tt/6 instead of 
(79) into the r.h.s. of (78), the latter reduces to (77) (with a = 
0). In the nonrelativistic limit both the time-symmetric and time- 
asymmetric (purely confinement) models lead to the same relation: 

M c - 2m « 3 (J^L) ' (82) 

which is known from the nonrelativistic linear confinement model. 
Classical Re or ere traiectories from the ereneral time-asvmmetric model 
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These purely classical results give us the base to consideration 
of semiclassical quantization of the model. By analogy with WKB 
approximation method we put 

S = h(£+£, £ = 0,1,... (83) 

for the quantized internal momentum, and 

/ k r dr = 27rh(n r + ±), n r = 0, 1, . . . (84) 

for radial excitations; the integral runs over the classical phase tra- 
jectory. 

In the case of purely confinement model we have: 

J dr^/f(r, M, S) = nh{n r + §). (85) 

Using the oscillator approximation!'! we expand the function 
f(r, M, S) (59) about the circular orbit to first nonvanishing orders 
in AM = M — M c and Ar = r — r c . The result is as follows: 

f(r, M, S) w a 2 (M, S) - b 2 (S)(Ar) 2 , (86) 

where 



dM 



AM = M jf t^AM, (87) 
(M r 2 + 8m 2 ) ' v 1 



2/m_ d 2 /(r,M,S) 



2<9r 2 



27/3 2 M 4 



(88) 



4(M c 2 -4m 2 )(M c 2 + 8m 2 )' 
and the function M C (S) is defined in (64). Then the integral in l.h.s. 
of (85) is easily calculated: 

a/b _ 2 

y d(Ar)v/a 2 -6 2 (Ar) 2 = |y. (89) 

—a/6 

Using (85)-(89) and assuming that AM is small compared to M c we 
obtain for M 2 the expression: 



{ ^/Zflh M 2 + 8m 2 1 
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which together with (83) and the definition (64) of the function 
M C (S) describes the leading (for n r = 0) and daughter (for n r = 
1,2,...) Regge trajectories. 

Similarly, in the case of general model we obtain the Regge tra- 
jectories determined in the implicit form by the equations 

M 2 = M C 2 (A) {1 + g$(A) (n r + h)} , (91) 



(I A ) = y [1 + 2A + ^(1 + A) 2 ] [3 + (l + ^)(l + A) 2 + 3Kl + A) 4 ] 

(l + A) 2 [l + i/(l + A)WA(2 + A) ' 1 } 



and (75), (83). At large t these trajectories reduce to linear ones, 

Ml « 6V3/3K{£ + n r + 1) + 6 (m 2 - 3a/3) , (93) 



so that the daughters are parallel to leading trajectory. Moreover, 
states of unit internal momentum differences form into degenerate 
towers at a given mass. This tower structure is of interest for the 
meson spectroscopy, as it is intimated in Ref. |6|. The number of rel- 
ativists potential models based on single-particle wave equations@S0 
as well as two-particle models with oscillator interaction0'ini lead to 
degeneracy of £ + 2n r type, but not of £ + n r type.^ The latter cannot 
be reproduced by single-particle relativistic models with the vector 
and scalar potentials, as it is shown in Ref. 

Figures 5 and 6 present two examples of semiclassical Regge tra- 
jectories which are characteristic for heavy and light mesons respec- 
tively. Trajectories in figures 5a and 6a are calculated in the oscillator 
approximation which is good for n r <C L Thus curved segments of 
daughters at n r > £ are not sure. This is evident by comparison 
to trajectories of the same case which are shown in figures 5b and 
6b. They are obtained by the numerical solving of the equation (84) 
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Figure 5: Semiclassical Regge trajectories; m = 1.27 GeV, a = 0.5, (3 = 0.2 GeV 2 ; 
a) oscillator approximation; b) numerical solution. 




Figure 6: Semiclassical Regge trajectories; m = 0.005 GeV, a = 0.8, f3 = 0.2 GeV 2 ; 
a) oscillator approximation; b) numerical solution. 
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where the integral in l.h.s. of (84) runs over phase trajectories of gen- 
eral model. We note that due to (83) phase trajectories in this case 
correspond to S > (see figure 1), and thus they are free of critical 
points discussed in Subsection 4.3. It is remarkable that numerical 
solutions for mass spectrum is well described by the asymptotic for- 
mula (93) even at small I. This is especially concerned with the case 
of light mesons (figure 6b). 

6. Conclusion 

In the present paper we have traced the relation between the nonstan- 
dard classical Yang-Mills field theory which arises from the consid- 
eration of QCD in the infrared regional and the classical relativistic 
two-particle models with confinement interaction formulated in the 
framework of Fokker formalism. It is notable that the use of Abelian 
potentials following from NT provides the confining interaction of 
particles. The time-symmetric (purely confinement) model turns out 
known in the literature^^ where, although, it was constructed as a 
priori action-at-a-distance model. 

The present time-asymmetric confinement model could be re- 
garded as a classical relativisation of the primitive quarkonium 
model. It has a number of features which are expected for mod- 
els of this kind but which usually are not realized together. 

1. The model is a self-consistent relativistic two-particle model. It 
allows the Lagrangian and Hamiltonian formulations. The quantities 
in terms of which the model is built have clear physical meaning. 
In the case S > solutions of this model lead to smooth timelike 
particle world lines. If S = 0, the collision critical points occur 

in ™rnir>n nort icloa roQr>n tlio cndprl nf lirrnf Tn f nia pqoc alflinnrrln 
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world line can be smoothly continued as well. In both cases a particle 
motion is bound. 

2. Estimates of Regge trajectory from classical mechanics shows 
that it has a proper asymptote while the corresponding nonrela- 
tivistic potential is linear. This feature is not derivable from non- 
relativistic models. The parameters of a linear rise following from 
the time-symmetric and time- asymmetric models differ from one an- 
other by near 5 %. One can hope that other long-range effects 
which should follow from the forthcoming study of purely retarded, 
time-symmetric and time-asymmetric confinement interactions differ 
slightly as well. 

The semiclassical consideration leads to the interesting degenerate 
tower structure of meson spectrum which probably exists in nature. 

3. The interpretation of an interaction in terms of some field the- 
ory is very important in RDIT. Hopefully, the knowledge of field 
equations and corresponding variational principle underlying the 
model allows to include properly into consideration spinning par- 
ticles and then to construct the quantum-mechanical description. 
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Appendix 

A. Green's functions for field equation of NT. 

Let us consider the equation (20) (here we omit unessential Greece 
indices and the factor k 2 ), 

U 2 A{x) = 4:7rJ(x). (A.l) 

It can be recast into the set of two d'Alembert equations, 

uA{x) = B{x), (A.2) 

UB(x) = 47rJ(x). (A.3) 
Solving them yields the following formal expression for A(x): 



A = —D v * B = — D v * Drf * J, (A.4) 



—D n * B = — 

where 77 and rf can independently take values +1, -1, or 0. 

Since the convolution of distributions is not guaranteed to be a 
well defined operation,^ we have to examine all possible combina- 
tions of 77 and rf . Actually, it is sufficient to consider the cases 
rj = ±1, 7/ = ±1; other cases where rj = or /and 7/ = reduce to 
the previous ones due to the linearity of equations (A.2), (A.3) and 
the equality Dq = \{D + + D ). 

Let us write down the expression (A.4) in the explicit form, 

A(x) = i- Jd 4 y Jd 4 zD v (y)D v ,(z)J(x -y-z). (A.5) 
Representing D v in the form 

D v (x) = 2e(f]x°)S(x 2 ) = S(x° - t/|x|)/|x|, 77 = ±1, (A.6) 
where x = {x\ i = 1, 2, 3}, yields for (A.5) the expression 
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or, in terms of new variables u = y + z, v = y — z, 
A(x) = -!- Jd 3 u Jd 3 v x 



J(x° — ^T]\u + V 




w\ 


u 


- V 


, x - u) 




u + v| 


|u - 


- V 





2/1 . ' 1 M 2/l " — ^. (A.8) 



Let us calculate the internal integral over d 3 v in r.h.s. of (A.8). 
Expressing Cartesian coordinates vi, i>2, v% of v in terms of ellipsoidal 
coordinates, 



Vi 



= |u|\/(cr 2 — 1)(1 — T 2 ) COS<^, 



v 2 = |u|v(cr 2 - 1)(1 - r 2 ) siny?, 

V3 = |u|crr, (A. 9) 

(J > 1 > r > -1, < p < 2tt, 



we obtain 

oo 1 2ir 

A(x) = — Jd 3 u | u| yiicr y<ir 



x 

1 -1 



xj(a»-„|u|x{*' J = ^),x-u). (A.10) 

If rf = —77 this integral diverges due the factor f^da. In the case 
rf = rj it reads: 

00 

A(aj) = i yii 3 ?i |u| JdaJ(x° - rj\u\a, x - u). (A.ll) 

1 

Using the change of the variable a — > w° = r/|u|cr we obtain the final 
expression for A(x): 

A(x) = \ Jd 4 ue(r]u )e(u 2 )J(x - u). (A.12) 

It follows from (A.12) that fundamental solutions to the equation 
(A.l) are: 
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Since their supports are the interiors of future- and past-oriented 
light cones, these distributions are the retarded and the advanced 
Green's functions of the equations (A.l). The time-symmetric 
Green's function is constructed by linearity, 

Eo(x) = \{E + {x) + E-(x)) = Mx 2 ). (A.14) 

Eqs. (A.13), (A.14) are unified in eq. (21). 

We note that some complex fundamental solution to (A.l) is ob- 
tained by means of another technique in Ref. 132, 



E c (x) = ^ ln(x 2 ± iO). (A.15) 

Its real part, 

m c (x) = \(Q(x 2 ) - 1), (A.16) 

coincides with the symmetric Green's function (A.15) up to a con- 
stant (which is the solution of homogeneous equation). This solu- 
tion can be considered as the analogue of the Feynman propagator 
in QED, 

D cW = i t\ - m ' ( A - 1? ) 

7T{X Z ± lO) 

real part of which is the symmetric Green's function of d'Alembert 
equation. 

B. Analysis of DPM for general model. 

Let us introduce the dimensionless positive quantities: 

2a M S a(3 

Mr 2m a m z 

and functions: 

/±«,M,*) = l + 2£- ^ 2 h ± (^fji,<T), (B.2) 
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£ ( i + a ± (f--j C+< * 



2^2 



(B.4) 
(B.5) 



which are related to (70)-(73) as follows: 

h±&fr&) = -^h ± {r,M,S) 
16 

4 

9-(£,V,v) = Jpd(r,M,S) 
Then DPM is determined by conditions: 



(B.6) 
(B.7) 
(B.8) 
(B.9) 



>0, (B.10) 

A = (l + O 2 -/-tt,M^)>0. (B.ll) 

Although ^ is positive by definition, it is useful to consider the func- 
tions (B.2)-(B.5) of £ for £ G H. 

First of all we consider the condition (B.10). It follows from (B.5) 
and (B.4) that 

g-(0>o i £e[(W/4 (B.12) 



d(Z) >?-{£), 



f G (-cx>,0)U(i//jU 2 ,oo), 
£ = 0,i//ju 2 , 
£g(0,^ 2 ). 



(B.13) 
(B.14) 
(B.15) 



Thus from (B.12)-(B.14) we obtain 

d(fl >0, £e (-cx>,0]U[i//ju 2 ,oo). 



(B.16) 



Using the equivalent form of d, 
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we obviously have 



d(0 > 

d(0 < 



£ G (—00, —1) U (0, v I fi 2 ). 
£ = -1,0,i//ju 2 , 

£g (-l,0)u(^ 2 ,oo), 



(B.18) 
(B.19) 
(B.20) 



and then, from (B.16) and (B.18). 



d(fl > 0, (ER. 



(B.21) 



Hence the functions h±(£), /±(£) (as well as h±(r), f±(r)) are real. 
Taking into account (B.3), (B.12)-(B.15), and (B.21) one obtains 



MO > 0, 
M0 = o, 
MO < 0, 



^G (O,^ 2 ), 
£ = <W// 2 , 

£G (-oo,0)U(i//jU 2 ,oo) 



(B.22) 
(B.23) 
(B.24) 



Thus, using (B.2) and (B.22)-(B.24) we conclude that 

/+(£)> 0, ^G [-l/2,oo). 



(B.25) 



We note that the function /+(£) is smooth at £ = i^//i 2 . 

In order to clarify the behaviour of the function /_(£) for £ > 
let us consider the function /+/_. It can be presented in the form: 



/ + (0/-(0 = ^ 1± ^n(0, 



where 

n(fl = (1 + 20 



(B.26) 



<j 2 (1 + 0£ 3 (B.27) 



is 4th-order polynomial (in terms of original quantities it is written 
down in l.h.s. of (74)). It is evident that 
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Moreover, it is easy to examine that 



n(fl < o, 
n(-i)>o. 



(-oo,-l] u[o,^ 2 ] 



(B.29) 
(B.30) 



Thus II(£) has two negative roots, £1 and £2, — 1 < £1 < — 1/2 < 
£2 < 0, which exist at arbitrary (positive) values of fi, a, and v. The 
number of positive roots, should they exist, is not more than two. 
Let us note that one can choose the sufficiently large value of a such 
that II(£) > 0, £ > 0. On the other hand, at S = there exists 
£ + > such that n(£) > 0, £ > £+. Thus, given ^ and i^, two other 
roots, £3 and £4, are positive for sufficiently small values of cr, and 
v/fj 2 < £3 < £4- In this case we have: 



Hereafter we restrict all functions on £ > 0. Using (B.26), (B.31), 
(B.32), and (B.25) one concludes that 



Now let us consider the condition (B.ll) which can be presented 
in the following form: 



n(fl > 0, 
n(fl < 0, 



£e (-00,6) u (6,6) u (£4,00). 



(B.31) 
(B.32) 




£e (0,i///* 2 )u 



(B.33) 
(B.34) 



A = . 2 , cr) + y/d{^, fl, a) > 0. 




(B.35) 



Using (B.18), (B.20) and the evident inequality 

p+(0 >0, £G W/« 2 ,oo), 



(B.36) 



one concludes that 



A > 



£ G (Wu 2 ,oc) 



(B.37) 
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Finally, taking into account (B.33), (B.34), (B.37), and (B.38), 
we find DPM: £3 < £ < £4, i.e., r\ < r < where T\ = 2a /(M £4), 
V2 = 2a/(M^). It disappear if £3 — > £4. The degenerated case 
£3 = £4 = £ c , where £ c satisfies the set of equations: 

n« c ) = 0, n'Kc) = 0, (B.39) 

corresponds to circular orbit motion. The set (B.39) in this case can 
be considered as a relation between M and S which is presented by 
(75), whereA = l/£ c . 
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